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We develop a kinetic theory of Brownian particles with long and short range interactions. Since 
the particles are in contact with a thermal bath fixing the temperature T, they are described by the 
canonical ensemble. We consider both overdamped and inertial models. In the overdamped limit, 
the evolution of the spatial density is governed by the generalized mean field Smoluchowski equation 
including a mean field potential due to long-range interactions and a generically nonlinear barotropic 
pressure due to short-range interactions. This equation describes various physical systems such as 
self-gravitating Brownian particles (Smoluchowski-Poisson system), bacterial populations experi- 
encing chemotaxis (Keller-Segel model) and colloidal particles with capillary interactions. We also 
take into account the inertia of the particles and derive corresponding kinetic and hydrodynamic 
equations generalizing the usual Kramers, Jeans, Euler and Cattaneo equations. For each model, we 
provide the corresponding form of free energy and establish the //-theorem and the virial theorem. 
Finally, we show that the same hydrodynamic equations are obtained in the context of nonlinear 
mean field Fokker-Planck equations associated with generalized thermodynamics. However, in that 
case, the nonlinear pressure is due to the bias in the transition probabilities from one state to 
the other leading to non-Boltzmannian distributions while in the former case the distribution is 
Boltzmannian but the nonlinear pressure arises from the two-body correlation function induced by 
the short-range potential of interaction. As a whole, our paper develops connections between the 
topics of long-range interactions, short-range interactions, nonlinear mean field Fokker-Planck equa- 
tions and generalized thermodynamics. It also justifies from a kinetic theory based on microscopic 
processes, the basic equations that were introduced phenomenologically in gravitational Brownian 
dynamics, chemotaxis and colloidal suspensions with attractive interactions. 



I. INTRODUCTION 

In the last ten years, the dynamics and thermodynamics of systems with long-range interactions (LRI) has been 
a subject of active research [iJtl. Systems with long-range interactions are numerous in nature and concern, for 
example, self-gravitating systems (galaxies and globular clusters), two-dimensional turbulence (vortices and jets), 
non- neutral plasmas, free electron lasers (FEL), and toy models such as the Hamiltonian Mean Field (HMF) model. 
In these systems, the interaction potential u(r) decays with a rate slower than l/r d at large distances, where d is the 
dimension of space (these potentials are sometimes called "non-integrable"). As a result, each particle interacts with 
far away particles (i.e. the interaction is not restricted to nearest neighbors) and the energy is non- additive. This 
can lead to striking properties (absent in systems with short-range interactions) such as inequivalence of statistical 
ensembles and negative specific heats in the microcanonical ensemble 

In a series of papers [Ill4l6| . we have developed a general kinetic theory of systems with long-range interactions. 
We have considered both isolated Hamiltonian systems described by the microcanonical ensemble and dissipativc 
Brownian systems described by the canonical ensemble. In the first case, the energy E is conserved while in the latter 
case, the system is in contact with a thermal bath fixing the temperature T . For systems with long-range interactions, 
the mean field approximation becomes exact in a proper thermodynamic limit N — > +00 where the coupling constant 
scales like k ~ 1/N while the volume of the system remains finite fl7Bl9l ]. For Hamiltonian systems, the dynamics of 
the one-body distribution function is described, in the N — > +00 limit, by the Vlasov equation [2fj |. This corresponds 
to a collisionless regime. A more general kinetic equation taking into account finite N effects can be obtained at the 
order 0{1/N) [la - [l5j . It describes the collisional regime. For spatially homogeneous systems, this kinetic equation 
corresponds to the standard Landau equation (when collective effects are neglected) |21| or to the Lenard-Balescu 
equation (when collective effects are taken into account) [H,Sll- For Brownian systems, the dynamics of the one-body 
distribution function is described, in the A~ — > +00 limit, by the mean field Kramers equation ff2| . In the strong 
friction limit £ — > +00, it reduces to the mean field Smoluchowski equation. In Refs. [HI Hj, [25| , we have developed a 
hydrodynamics of Brownian particles in interaction. The first two moments of the hierarchy of hydrodynamic equations 
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are the damped Jeans equations. If we implement a Local Thermodynamic Equilibrium (LTE) approximation, we 
obtain the damped Euler equations. The validity of this approximation has been investigated in |25| in the case of 
simple models. If we neglect the velocity tensor in the damped Euler equations, we get a Cattaneo-type equation for 
the density taking into account memory effects. Finally, in Ref. (l6j | we have developed a theory of fluctuations in 
which the previous deterministic partial differential equations are replaced by stochastic partial differential equations 
including a multiplicative noise term depending on position and time. Fluctuations are important (i) when the number 
of particles is small, and (ii) close to a critical point. When the system possesses metastable states, the fluctuations 
can induce random transitions from one state to the other. These kinetic and hydrodynamic equations, involving 
long-range interactions, can model various physical systems such as self-gravitating Brownian particles [26| . bacterial 
populations experiencing chemotaxis j27j and colloids at a fluid interface driven by attractive capillary interactions [28| . 
These models were initially studied in the strong friction limit in which the dynamics of t he p articles is overdamped 
but, later, it was realized (e.g. in chemotaxis) that inertial effects [29l - [3l| and fluctuations |32| can play an important 
role so that hydrodynamic and stochastic models have also been introduced. 

The d yna mical evolution of Brownian particles in interaction is also studied in the physics of simple liquids and 
colloids [331 ] . In that case, the interactions are short-ranged and the mean-field approximation is not valid. The 
usual approach to take correlation functions into account is based on the Density Functional Theory (DFT) [34| 
and on the Dynamical Density Functional Theory (DDFT) [35|. The correlations between particles induced by 
a short-range potential of interaction lead to an excess pressure p ex (r,t) with respect to the ideal gas pressure 
Pid(r,t) = p(r,t)kBT/m. Deterministic and stochastic models of interacting Brownian particles have been developed 
in that context . 

For certain systems, it is important to take into account both long-range and short-range interactions. For example, 
due to the attractive long-range interaction, self-gravitating Brownian particles, bacterial populations and colloids 
driven by attractive capillary interactions can collapse. In that case, the central part of the system becomes very dense. 
In the absence of short-range interactions, the collapse generically leads to the formation of Dirac peaks (26| . However, 
these peaks are unphysical and, in practice, the density profile is regularized by small-scale constraints. These small- 
scale constraints can be due to finite size effects (the particles always have a finite size and cannot interpenetrate), 
steric hindrance, short-range interactions and, ultimately, quantum mechanics (Pauli exclusion principle). These 
interactions will come into play when the system becomes dense enough. Their effect is to provide a nonlinear 
pressure that will halt the collapse and lead to a well-defined equilibrium state. An example of this rcgularization is 
provided by a gas of self-gravitating fermions in which gravitational collapse is balanced by the pressure force arising 
from the Pauli exclusion principle [1(1 H3 ■ 

The aim of this paper is to present a unified kinetic theory of Brownian systems that takes into account both long 
and short range interactions. This will provide a precise justification, from a microscopic theory, of the kinetic and 
hydrodynamic models that have been introduced phenomenologically to describe self-gravitating Brownian particles 
[26| . chemotaxis 27 1 and colloids with capillary interactions (28j. This will also make the bridge between the physics of 
long-range interactions Q and the physics of simple liquids with short-range interactions (33l . 34 1 . The approximations 



made in these two topics are radically different. For long-range interacting (LRI) systems, the mean field approxi- 
mation applies 0] while for short-range interacting (SRI) systems, the correlations between particles are crucial and 
must be taken into account f33l H3] • We shall also develop some connections with nonlinear mean field Fokker-Planck 
equations (see reviews [HI, l39j ) based on generalized thermodynamics pol - li^ . In the context of generalized ther- 
modynamics, nonlinear Fokker-Planck equations arise when the transition probabilities from one state to the other 
depend in a non trivial manner on the occupancy of the starting and arrival states [43j . This kinetical interaction 
principle (KIP) takes into account microscopic constraints that affect the dynamics of the particles at small scales. 
Interestingly, this yields the same type of hydrodynamic equations [39l . [44| as in the case of short-range interactions 
although the justification of the nonlinear pressure is different. In the case of generalized thermodynamics, the non- 
linear pressure is due to the bias in the transition probabilities that leads to non-Boltzmannian distributions while in 
the DFT and DDFT used in the physics of liquids the distribution is Boltzmannian but the nonlinear pressure arises 
from the two-body correlation function induced by the short-range potential of interaction. These two approaches 
therefore take into account microscopic constraints in a different manner. However, it is interesting to find some 
connections between the hydrodynamic (macroscopic) equations although the kinetic equations are different. This 
shows that these approaches are complementary. 

As a whole, our paper develops connections between the topics of long-range interactions, short-range interactions, 
nonlinear mean field Fokker-Planck equations and generalized thermodynamics. The applications of our kinetic theory 
concern various physical systems such as self-gravitating Brownian particles, chemotaxis of bacterial populations, 
colloidal particles with capillary interactions, and probably others. Since our paper brings together several topics, 
and is addressed to an audience with different backgrounds, it is necessary to briefly review the most important results 
of each topic to make the paper self-contained. However, the kinetic equations that we obtain are new and generalize 
those obtained separately in the physics of long-range interactions and in the physics of liquids. 
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II. STATISTICAL EQUILIBRIUM STATE OF BROWNIAN SYSTEMS 
A. The Gibbs canonical equilibrium 

We consider a system of N Brownian particles with identical mass m interacting via a potential U(rx, r^r). We 
assume that the potential is of the form 

U(ti, ...,r N ) = m 2 ^2u LR (\r l - r 3 |) + m 2 ^usR(|rj - r,-|) + $ e2;t (r i ), (1) 

i<j i<j i 

where ulr is a long-range binary potential, usr is a short-range binary potential and $ ex t is an external potential. 
The Hamiltonian is 

- 1 

H = J2^iv 2 + U(r u ...,r N ), (2) 

i=l 

where the first term is the kinetic energy K and the second the potential energy U. These Brownian particles are in 
contact with a heat bath with temperature T so that they are described by the canonical ensemble. The statistical 
equilibrium state is given by the Gibbs canonical distribution 

P n (t u v x ,...,t n ,v n ) = -l— e -P^ m ^+ u ^-"^\ (3) 

^tot(P) 

where j3 = l/(fcsT) is the inverse temperature. The iV-body distribution Pjv(ri, vi, rjv, vjv) gives the 
probability density that the first particle is in (ri, Vi), the second in (r2, V2) etc. The normaliza- 
tion condition J P^dr 1 dv 1 ...dr^dvN = 1 leads to the expression of the partition function: Z tot (f3) = 
J exp[— /?H (ri, Vi..., rjv, Vjv)]rfrirfvi...drjvrfvjv- We introduce the free energy functional 

Ftot [Pn] = E tot [P N ] - TStot [Pn] , (4) 

where 



Stot = -ks J Pn In Pn dr 1 dw x ...dr N dv Nl (5) 

is the entropy and 

(H) = J PnH dridvi...dr^dvN, (6) 



E,, 



is the average energy. The canonical ./V-body distribution ([3]) minimizes the free energy Ftot[PN] at fixed normalization. 
Furthermore, the value of the free energy at equilibrium, obtained by substituting the Gibbs distribution §S§ in Eqs. 
(O and ©, is F to t{P) = — (V/?) hiZ tot (/3). The average energy at equilibrium is E tot (f3) = — dm Z tot /df3. The 
fluctuations of energy are ((AH) 2 ) = C/(k B /3 2 ) where C = 8E tot /dT = —kn(3 2 dEt„t/d0 is the specific heat. This 
relation shows that the specific heat is always positive in the canonical ensemble [45| . 

From Eq. we see that the velocity dependence of the TV-body distribution is Gaussian. Therefore, the average 
kinetic energy is (K) = dNksT/2 (where d is the dimension of space) just like in a non-interacting gas. In the 
following, we shall mainly focus on the configurational part of the distribution function 

P N (r 1 ,...,r N ) = ^e-^-^\ (7) 

which contains the non-trivial information on the system. The normalization condition J Pjv^ r i---rfrjv = 1 leads to 
the expression of the partition function: Z(f3) = J exp[— pjJ(r\, rjv)]rfri...drjv- Due to the Gaussian nature of the 
velocity distribution, we have 

P N (r u v 1 ,...,r N ,v N )=(^\ e-^-^P^n,..,!^). (8) 

Comparing Eqs. ©, © and ©, we find that Z tot (P) = (2tt/ f3m) dN / 2 Z((3). We introduce the configurational free 
energy 

F[P N ] = E[P N ]-TS[P N ], (9) 
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where 



is the configurational entropy and 



S = -k B J P N lnP N dr 1 ...dr Nl (10) 

E = (U) = J P N U dri...dr N , (11) 

is the average potential energy. The canonical ./V-body distribution (JTJ) minimizes F[Pn] at fixed normalization. 
Furthermore, the value of the free energy at equilibrium, obtained by substituting Eq. (JT)) in Eqs. (|10[) and (|11[) . is 
F(f3) = — (1//3) hxZ(/3). The average energy and the fluctuations of energy are given by expressions similar to those 
given above. 

Remark: the free energy F[Pn] does not always have a minimum. This is the case in particular for the gravitational 
interaction in d = 3 due to the phenomenon of gravitational collapse. The strict statistical equilibrium state of a 
self-gravitating gas in the canonical ensemble is a Dirac peak containing all the particles (46l - |48j . Such a configuration 
makes the free energy diverge to — oo due to the (algebraic) divergence of the potential energy that cannot be 
compensated by the (logarithmic) divergence of the entropy in the other direction (47| . However, there can exist 
metastable states in the form of gaseous configurations that have very long lifetimes, scaling like e N [49| . These 
mctastable states are local minima of the mean field free energy functional F[f] defined by Eq. (p3|) . See Appendix 
|B]for some comments about the importance of metastable states. 



B. The Yvon-Born-Green (YBG) hierarchy 

We introduce the reduced probability distributions 

Pj(ri,...,Tj) = J P/v(ri, ...,r N ) dr j+1 ...dr N . (12) 

Differentiating the defining relation (fT2"j) for Pj and using Eq. ([7]) , we obtain the YBG hierarchy of equations [ill, HH : 

9Pj , x „ o „ / x du\ i 

^i(r lj ...,r i ) = -/3m 2 P i (r lj ...,r J )X;-^' 



i=2 



-Pm 2 (N-j) J Pj+1 { ri ,,.., rj+l )^^dv j+1 -PmP j {v 1 ,...,v j )^, (13) 

where u = Ulr + usr denotes the total binary potential of interaction and we have noted Uij for u(|rj — Yj\). The 
first equation of the hierarchy is 

^(r 1 ) = -/3m 2 (7V-l) / P 2 {j uV2 ) d -^idv 2 - PmP.iv,)^^. (14) 
OTi J or i or i 

If we introduce the local density p(r) = NmPi(r) and the two-body distribution function p2(r,r') = N(N — 
l)m 2 P 2 (r, r'), the first equation of the YBG hierarchy becomes 

Ml V p(r) = - / p 2 (r,r / )Vu(|r-r'|)dr'-p(r)V* est (r), (15) 
m J 

where we used the fact that the particles are identical. This equation determines the equilibrium density profile 
p(r) when the two-body correlation function is known. Since the velocity distribution is Gaussian, the equilibrium 
distribution function /(r,v) = NmPi(r,v) is given by 



/3 m 



d/2 



On the other hand, the free energy functionals ((U) and © defined in Sec. Ill Al can be written 

Ftot[P N ] = \j fv 2 drdv+^ J p 2 (r,r')u(\r-r'\)drdr' + J p<P ext dr + k B T J P ff \n.P N dr 1 dv 1 ...dr N dv N , (17) 



and 



F[Pn] = \J p 2 (r,r>(|r-r'|)drdr' + j p<P ext dr + k B T J P N lnP I fdr 1 ...dr N . (18) 
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C. The virial theorem 

Taking the scalar product of Eq. (|1 5|) with r, integrating over the entire domain and integrating by parts (assuming 
that the boundary terms can be neglected), we obtain the exact virial theorem 

dNk B T + V LR + V SR + V ext =0, (19) 

where 

V L r = - J p 2 (r, r')r • Vu LR (v, r') dvdv\ (20) 

V SR = - J p 2 (ry)r-VusR(r,r')dTdv', (21) 

V ext = - J p(r)r ■ V$> ext (r) dr. (22) 

The first term is twice the average kinetic energy (A'). The second term is the virial of the long-range interaction, 
the third term is the virial of the short-range interaction and the fourth term is the virial of the external force. 

D. Long-range interactions: mean field approximation 

We first consider a purely long-range interaction. In that case, it has been established rigorously [l8| that, in a 
proper thermodynamic limit TV — > +oo [74| , the mean field approximation is exact: the TV-body distribution function 
is a product of TV one-body distribution functions 

P N ( ri ,...,r N ) = Pi(r 1 )...P 1 (r JV ). (23) 

More precisely, it can be shown that the non trivial correlation functions P' n of order n scale like TV~(™ -1 ) [ll|. In 
particular, P2(i"i,r 2 ) = Pi(ri)Pi(r 2 ) + 0(1/TV). Therefore, we can make an expansion of the equations of the YBG 
hierarchy in powers of the small parameter 1/TV [ll[. Here, we limit ourselves to the limit TV — > +oo so that 

p 2 (v,v') = p(r)p(r'). (24) 

In the mean field approximation, the first equation (fT5j) of the YBG hierarchy becomes 

T Vp(r) = -p(r)V / p(r')u(\r - r'\)dr'. (25) 



m 

After integration, this can be written in the form of a mean field Boltzmann distribution 

p{r)^Ae-^ {r \ (26) 

where 

$(r) = J p{v')u LR {\v-T'\)dv', (27) 

is the mean field potential produced self-consistently by the smooth distribution of particles. Therefore, the equilibrium 
density profile of the particles is determined by an integrodifferential equation (f25|) . Since the velocity distribution is 
Gaussian, the distribution function in phase space is given by the mean-field Maxwcll-Boltzmann distribution 

/(r,v) = A'e-^^+^M]. (28) 
Equation (|25[) can be written as a condition of mean field hydrostatic equilibrium 

Vf>(r) = -p(r)V$(r), (29) 
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with an equation of state 



p(r) = p(r)%^. (30) 



Therefore, in the mean held approximation, the local equation of state of a Brownian gas coincides with the isothermal 
equation of state. Note that the pressure appearing in Eq. (j2l)j) is the kinetic pressure defined by 

PfeinW = 3 [ fv 2 dv. (31) 



The virial theorem is given by Eq. (fT9]) where Vsr = V ext = and where Vlr is replaced by the mean field virial 

V = - ( pr- V$dr. (32) 



For the gravitational potential in d dimensions, which is solution of the Poisson equation A$ = SdGp, the mean field 
virial is V = (d — 2)W for d ^ 2 (where W = ^ f p$dr is the mean field potential energy) and V = —GM 2 /2 for 
d = 2 [24[ . Substituting Eq. (J23J) in Eq. (|18[) , and introducing the mean field potential (|27[) , we obtain the free energy 
functional 

F[p] = E[p] - TS[p] = - I p$>dv + k B T [ In -?-dr, (33) 
2 J J m to 

up to an additive constant term —NksT In iV. Similarly, in the mean field approximation, the free energy (|17p reduces 
to the form 



F[f] = E[f] - TS[f] = \ ( fv 2 drdv + \ [ p4>dr + k B T f — In ^-drdv, 
2 J 2 J J m m 



(34) 



up to an additive constant term — iVfcsTlniV. It can be shown rigorously that the equilibrium density minimizes F[p] 
at fixed mass [l8j . This yields the mean field Boltzmann distribution (|26[) . Similarly, the equilibrium distribution 
function minimizes F[f] at fixed mass. This yields the mean field Maxwell-Boltzmann distribution (|28j) . 

Remark: It has to be noted that the mean field equation (|25p may have several solutions for a given value of mass 
M and temperature T (see, e.g. |47|], for self-gravitating systems). Only (local) minima must be selected. The global 
minimum corresponds to the strict equilibrium state. A local (but not global) minimum corresponds to a metastablc 
state. Saddle points must be rejected because they are "unstable" for some perturbations. 

E. Long and short-range interactions: Density Functional Theory (DFT) 

1. Exact results 

We now consider a system of Brownian particles at temperature T with long-range and short-range interactions (for 
the sake of generality we also assume that they evolve in a fixed external potential). A central result in the theory of 
fluids [33l [34l | is that, even if there exists non trivial correlations between the particles, the equilibrium density profile 
p(r) minimizes a free energy F[p] at fixed mass. This free energy can be written as 

F[p] = [ p$ ext dr + 1 / p<P dr + k B T [ ^ In ^dr + F ex [p\. (35) 
J 2 J J to m 

The first term is the potential energy associated with the external potential. The second term is the mean field 
potential energy associated with the long-range potential of interaction (see Sec. Ill D|) . The third term is the free 
energy of the ideal gas. Finally, the fourth term is the excess free energy .F ear [/o]. This is a non-trivial functional 
determined by the short-range interactions. All the difficulty in the theory of fluids is to find some approximate forms 
of this functional. Once this functional is known, the density profile, as well as all the 71-point correlation functions, 
can be obtained via functional differentiation. Inversely, the excess free energy is often obtained from the study of 
the correlation functions. The excess free energy F ex is known exactly only in a few particular cases, but very good 
approximations can be devised in more general cases [33|. 
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The fact that the density profile minimizes a free energy functional at fixed mass implies that its first constrained 
variations vanish. Writing SF + aTSM = 0, where a is a Lagrange multiplier taking into account the conservation of 
mass, we get 



Sp 

Taking the gradient of this expression, we obtain 



aT = 0. (36) 



Sp 

With the decomposition (f3"5"j) . these equations can be written 

k B T „, , 



v(? S )=o. 



[1 + ln(p(r)/m)] + $(r) + $ ea!t (r) + + aT = 0, (38) 

dp 



and 



k T 5F 

— Vp + pV — — +pV$ + pV$ ext = 0. (39) 
m op 

The equilibrium density profile is given by 

p(r) =Ae-M*+*" t+ ^). (40) 

We stress, however, that the r.h.s. depends on p(r) itself, so that Eq. P0|) is an integrodiffcrcntial equation. 

On the other hand, the first equation (|15|) of the YBG hierarchy including long-range and short-range interactions 

is 



k B T 



m 

Comparing Eqs. ([39]) and (|4Tj) , we obtain 



Vp = - J p 2 (r, v')Vu SR {r, r') dr> - pV$ - pV$ ext . (41) 



f SF 
J P2 (r,r')V US fl(r,r')dr' = p(r)V^[p(r) 



(42) 



This relation is exact at statistical equilibrium and is a central result in the theory of fluids j33l |34| . It relates the two- 
body correlation function to the excess free energy functional. Then, Eq. (jH]) can be viewed as an integrodiffcrcntial 
equation determining the equilibrium density profile once the excess free energy is known. 

2. Virial theorem 

Let us introduce the DFT virial 

V DFT = - J pr ■ V ^ dr. (43) 

At equilibrium, according to the exact identity (|42p. we have 

V SR = V DFT - (44) 
On the other hand, for the long-range interaction, we can make the mean field approximation and use 

V LR = V. (45) 
Therefore, the virial theorem is given by Eq. (fl~9| with Eqs. (|44[) and (|45l) . 
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3. Barotropic pressure 

In a fluid, the local pressure is of the form p = p(p,T). Since the temperature T is fixed in the case of Brownian 
particles (canonical description), the pressure is barotropic and we shall simply write p = p(p). In principle, the excess 
free ener gy F ex [p] can depend on the gradients of the density. This is particularly important for a fluid close to an 
interface [34|]. Here, we shall assume that the density varies on a distance that is large with respect to the range of 
intermolecular forces. This is the case if the density distribution is mainly due to long-range interactions, as we shall 
assume in the following. With this assumption, the free energy is of the form (see Appendix |A|) : 

F[p] = \jp<S>dv + J p$ ext dr + J pf dr - (46) 

The excess free energy is therefore 

Fex [p}= I P I" ^ir d Pi dv - k sT I - In -dr. 
J J PI J m m 

We note the relation 



(47) 



k T S F 

Vp(p) = — Vp + pV— — = Vp ld + V Pex , (48) 
m op 

where Pid(r) = p(r)kFjT/m is the ideal pressure law and p ex is the excess pressure due to short-range interactions. 
This relation can be used to determine the equation of state p{p) corresponding to the excess free energy F ex [p] and 
vice versa. For an ideal fluid (F ex = 0), we recover the perfect gas law p(r) = p(r)kBT/m. On the other hand, using 
the exact identity (|42|) . we find that 



Vp=-^-Vp(r)+ / p 2 (r,r')Vu5fl(r,r')dr'. 
m J 



(49) 



This equation relates the local pressure p(r), hence the equation of state p(p), to the two-body correlation function. 
Taking the scalar product of this relation with r, integrating over the entire domain and integrating by parts (assuming 
that boundary terms can be neglected), we obtain 

J pdr = Nk B T-^J p 2 (r,r')r- \7u SR {r,r r ) drdr' . (50) 
For a spatially uniform fluid, we recover the virial equation (33[: 

V = 1 ~ 9(0u'sR(m d dt, (51) 

where n = p/m = N/V is the number density and g(£) the radial correlation function defined by p2(r, r') = n 2 g(\r — 
r'|). 

4- Hydrostatic equilibrium 

For a free energy of the form (|4^]l , Eq. reduces to the condition of hydrostatic equilibrium 

Vp + pV$ + pV$ext = 0. (52) 
Comparing this relation with the first equation of the YBG hierarchy 

^Vp(r) + J p 2 (r,r')Vu SR (r,r')dr' + J p 2 (r, r')Vu LR (r, r') dr' + p(r)V$ e , t (r) = 0, (53) 

we note that the long-range interactions create a mean field force — pV^ while the short-range interactions create 
an excess pressure p e x(r) with respect to the ideal pressure law Pid(r) = p(r)fcsT/m. Long range and short range 
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interactions have therefore a very different influence on the system. On the other hand, for a free energy of the form 
flU, the DFT virial (g2]) takes the form 



V DFT =d j pdr - dNk B T. (54) 
Using Eqs. (|44|) . (j45]l and (f54|) . the virial theorem ([19]) can be written explicitly 

d / pdr + V + V ext = 0. (55) 



5. Weakly inhomogeneous systems 

In Sec. HIE 31 we have given a first simplified expression of the excess free energy. Here, we shall briefly mention 
another simplified expression that has been extensively studied in the physics of liquids (see, e.g. [HI, HH ) . If the 
density distribution varies slowly and exhibits small departures relative to some reference density p, we can expand 
the functional i^^fp] to second order in |p(r,r) — p| <C ~p, thereby obtaining 

F e A P ] = -\k B T J{p-p){v,t)c{\T-v'\^){ P -p){v',t)dvdv', (56) 

where c(r,p) denotes the Ornstcin-Zcrnikc direct correlation function in the homogeneous reference system 

There are several methods in the physics of liquids to obtain useful approximations of the direct correlation function 
[IH, hence of the functional (j51))) . Interestingly, we note that the functional ([55)1 has the same form (up to a shift in 
density [75| ) as the mean field free energy functional (f3"3")l with Eq. ([27)1 provided that we view the direct correlation 
function as an effective binary potential 

c(r) = -f3u eff (r). (58) 

This makes possible to apply the results obtained for mean field potentials to this particular situation by using the 
correspondence 



III. KINETIC THEORY OF BROWNIAN PARTICLES IN THE OVERDAMPED LIMIT 

A. BBGKY-like hierarchy 

In the overdamped limit, the dynamics of N Brownian particles in interaction is governed by the coupled stochastic 
equations [ijj: 

£— ^= V l U(r 1 ,...,r N ) + V2DR l {t), (59) 

at m 

where Ri(i) is a Gaussian white noise such that (Ri(t)) = and (Rf(t)R?(t')) = 5 i:j 5 a p5(t - t'). Here, i = 1, N 
label the particles and a = 1, ...,d the coordinates of space. These equations can be obtained from Eqs. ([75|) -(|76 ]) in 
the strong friction limit £ — > +oo or, equivalently, for large times t 3> In these limits, it is possible to neglect 
the inertial term in Eq. ([75]) leading directly to Eq. (|5"9"|) . The diffusion coefficient in phase space D is related to the 
friction coefficient £ and to the temperature T by the Einstein relation [53| : 

D=&**. (60) 
m 

In terms of the mobility fi = l/£m and of the diffusion coefficient in physical space = -D/£ 2 , the Einstein relation 
takes the form D* = fiksT. The time evolution of the iV-body distribution Pjv(ri, r^r, t) is governed by the iV-body 
Fokkcr-Planck equation 



A' 



dP N _ d 

1=1 



k B TdP N , i d , 

—Pn-z-U{ti, ...,rjv) 



(61) 
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This particular Fokker-Planck equation is called the TV-body Smoluchowski equation. The iV-body Smoluchowski 
equation monotonically decreases the free energy Indeed, a direct calculation yields the canonical .ff -theorem: 



N 

* = -E 

i=l 



1 



k B TdP d 



N 



£toPjv V. m d r i 



1 dU N 

— Pat— ) drx...dTN < 0. 



(62) 



For a steady state, F = 0, the term in parenthesis in Eq. (jf?2"j) . which is the diffusion current in the Smoluchowski 
equation (foTj) . vanishes. This leads to the Gibbs canonical distribution (JT]). Therefore, the Brownian gas described 
by the stochastic equations (|59p automatically relaxes towards the Gibbs distribution (if it is normalizable). 

It is easy to derive from Eq. (|6T|) the equivalent of the BBGKY hierarchy for the reduced distribution functions 
(1121). It reads (H: 



dP 

? dt 



E 

i=l 



9ri 



m 9ri 



^ E 



+ (JV - j)m / P 



9r, 



dr i+ i + Pj 



9r, 



(63) 



The stationary solutions of these equations coincide with the equations (fT3j) of the YBG hierarchy. The first equation 
of the BBGKY-like hierarchy is 



SPi d 



dt 



dri 



+ (N - l)m / P 2 — — dr 2 + Pi 



m dri 



dri 



dri 



(64) 



Introducing the local density p{r,t) = NmPi(r,t) and the two-body distribution function p 2 (r, r',t) = N(N 
l)m 2 P2(r, r', t), we obtain the exact Smoluchowski equation 



^dt 



k R T. 



Vp+ / p 2 (r,r',i)Vu(|r-r'|)rfr' + pV$ e 



(65) 



where we have used the fact that the particles are identical. The steady state of this equation returns the first YBG 
equation (fT5|) . 



B. Long-range interactions: mean field Smoluchowski equation 



For systems with long-range interactions, the mean field approximation (|23|) . extended out-of-cquilibrium, is exact 
when N — > +oo. In particular, we have 



p 2 {r,r',t) = p{r,t)p{r',t). 



(66) 



Substituting this relation in the first equation (|65|) of the BBGKY-like hierarchy, we obtain the mean field Smolu- 
chowski equation 



e |_v.(i£W.,, 



with 



$(r,t) = J p(r\t)u LR (\r-r'\)dr'. 
The free energy associated with this equation is given by Eq. f)33l) . 



(67) 



(68) 



C. Long and short-range interactions: Dynamical Density Functional Theory (DDFT) 



In the presence of long and short range interactions, the first equation (jf?5j) of the BBGKY-like hierarchy can be 
written 



e3e = v 



k B T 



m 



Vp+y p 2 (r,r',t)Vus R (\r-r'\)dr' + pV$ + pV<f> e 



(69) 
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where we have used the mean field approximation to evaluate the long-range interaction term. This equation is not 
closed since it depends on the two-body distribution function. In the dynamic density functional theory of fluids, the 
hierarchy is closed by making the approximation J3f|: 



J P2 (r,r' ,t)Vu SR (\r - r'\) dr' = p(r,t)V 



SF ec 
Sp 



■\p(r,t)], 



(70) 



where -/^[p] is the excess free energy calculated at equilibrium. This relation is exact at equilibrium (see Sec. Ill E|) 
and the approximation consists in extending it out-of-equilibrium with the actual density p(r, t) calculated at each 
time. This closure is equivalent to assuming that the two-body dynamic correlations are the same as those in an 
equilibrium fluid with the same one body density profile. Although it is not possible to ascertain the validity of this 
approximation in the general case, it has been observed for the systems considered that this approximation gives 
remarkable agreement with direct Brownian iV-body simulations. With this approximation, Eq. (|69[) becomes 



ksT _ „SF ex 
m dp 



(71) 



which is closed. If we consider a free energy of the form (|46j) , then using identity (|48[). the foregoing equation can be 
rewritten 



^ = V • (Vp + pV$ + pV<S> ext ) . 



(72) 



This is a generalized mean field Smoluchowski equation including a generically nonlinear barotropic pressure p(p) due 
to short-range interactions and a mean field potential $ = ulr * p (where * denotes the product of convolution) due 
to long-range interactions. This equation, coupled with a potential of the form A$ — fc 2 <f> = SdGp (where A: -1 is a 
screening length) arises in several physical problems such as self-gravitating Brownian particles |2(| , chemotaxis [27| , 
colloids with capillary interactions [28[, etc. By combining results issued from the physics of systems with long-range 
interactions [|| with those issued from the dynamic density functional theory of fluids (35|, we have here justified this 
equation from a microscopic model. 

Introducing the free energy functional (|35[) . we can write the generalized mean field Smoluchowski equation (1711) in 
the form 



dp 

at 



1 SF 



(73) 



This equation monotonically decreases the free energy functional ([55)) which plays therefore the role of a Lyapunov 
functional. Indeed, a straightforward calculation leads to the iJ-theorem appropriate to the canonical ensemble 



F 



SF dp 
~Sp~di' 



dr 



5F, 
Sp 



1 ^SF 

£ 5 p 



dr 



V^Vdr<0. 
dp J 



(74) 



For a steady state, F = 0, the last term in parenthesis vanishes so that SF/Sp is uniform. This leads to Eq. (|37l) . 
Therefore, a density profile p(r) is a steady state of the generalized mean field Smoluchowski equation ([73"]) iff it is a 
critical point of F at fixed mass. Furthermore, it can be shown that a steady state is linearly dynamically stable with 
respect to the generalized Smoluchowski equation ([73"]) iff it is a (local) minimum of F at fixed mass 38 , [3i| . This is 
consistent with the condition of thermodynamical equilibrium (see Sec. Ill E[) . If F is bounded from below (76[, we 
know from Lyapunov's direct method that the system will converge towards a (local) minimum of F at fixed mass M 
for t —> +oo. If several (local) minima exist (metastable states), the choice of the selected equilibrium will depend on 
a complicated notion of basin of attraction. 

Remark: Eq. (|73[) can be justified in a phenomcnological manner from the linear thermodynamics of Onsager if we 



interpret it as a continuity equation dtp + V • J = with a current J = — (l/£)V|p proportional to the gradient of a 
potential p(r) = SF/Sp that is uniform at equilibrium (see Eq. (1361 )). This is precisely the way in which this equation 
was initially introduced in the physics of liquids 0, Hj]; see also |38l. l39l| in a more general context. 
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IV. KINETIC THEORY OF BROWNIAN SYSTEMS TAKING INTO ACCOUNT INERTIAL EFFECTS 

A. BBGKY-like hierarchy 

We now take into account inertial effects and consider N Brownian particles in interaction described by the coupled 
stochastic equations [l2T |: 

It = (75) 

^ = -£v< - -ViU(v u ...,v N ) + V2DHi(t), (76) 
at m 

where — is a friction force, D the diffusion coefficient in phase space and Rj(t) a Gaussian white noise. The diffusion 
coefficient and the friction force are related to each other by the Einstein formula (f^Uj) . The system is described by 
the canonical ensemble where the temperature T measures the str eng th of the stochastic force (since D ~ T). The 
stochastic process (|75l) - (|76l) extends the classical Brownian model [531 ] to the case of particles in interaction. In this 
context, the friction is due to the presence of an inert gas and the stochastic force is due to classical Brownian motion, 
turbulence or any other stochastic effect. The friction and the noise can also mimic the overall influence of an external 
medium (not represented) with which the particles interact. This is the notion of "thermal bath" . The evolution of 
the N-body distribution function Pjv(ri,vi, rjy, vjv, t) is governed by the iV-body Fokker-Planck equation 

d P N dP N 8P N \ 9 ( k BTdP N \ 

where Fj = — — VjJ7(ri, ...,rjv) is the force by unit of mass acting on the z-th particle. This particular Fokker-Planck 
equation is called the N-hody Kramers equation. For £ = D = it reduces to the Liouville equation which governs 
the evolution of an isolated Hamiltonian system in the microcanonical ensemble [l2| • The TV-body Kramers equation 
decreases the free energy Q. Indeed, a direct calculation yields the canonical if -theorem: 

F=-V] [ lr( — ^ + Psvt) Vdvi-dr^ < 0. (78) 

For a steady state, F = 0, the term in parenthesis in Eq. (|78[) , which is the diffusion current in the Kramers equation 
(fTTj) . vanishes. Since d/dt = 0, the advective term (l.h.s.) in Eq. (|77|l must also vanish, independently. From these 
two requirements, we find that the stationary solution of the iV-body Fokker-Planck equation is the Gibbs canonical 
distribution ([3]). Therefore, the Brownian gas described by the stochastic equations (|75|) - (|76)) automatically relaxes 
towards the Gibbs distribution. 

It is easy to derive from Eq. ([771) the equivalent of the BBGKY hierarchy for the reduced distribution functions 
Pj(ri, vi, Yj, Vj, t). It reads fl2j : 

<>l'j . \ - dPj v - \ - dPj . x - /' On,,, . i dP j+1 

i=l i=l k=l,k=£i i=l 

B^ext 8Pa X d fk B TdP, 



. (79) 



2=1 2=1 X 

In particular, the first equation of the hierarchy is 

dPx , dPi d f A 9U ^J A 9Pl 

+ Vl ' n ( N ~ 1 ) m ^ / p 2(ri,vi,r 2 ,v 2 ,t)— — rfr 2 dv 2 — 

at oti avi J or i or i ovi 

d fk B TdP 1 \ , 

Introducing the distribution function f(r,v,t) = NmPi(r,v,t) and the two-body distribution function 
/ 2 (r, v, r', v') = iV(./V — l)m 2 P 2 (r, v, r', v'), we obtain the exact Kramers equation 



where we have used the fact that the particles are identical. 
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B. Long-range interactions: mean field Kramers equation 

For systems with long-range interactions, the mean field approximation is exact when N — > +00. In particular, we 
have 

/ 2 (r,v,r',v',t) = /(r,v,t)/(r',v',t). (82) 

Substituting this relation in the first equation (f8Tj) of the BBGKY-like hierarchy, we obtain the mean field Kramers 
equation 

^ +v .^Z_ V <I>-^-=£— • ( k -*L d J- +fv) (83) 
dt dr dv dv \ m dv J ' 

where $(r, t) is given by Eq. (|68p . For £ = D = 0, we recover the Vlasov equation which describes Hamiltonian 
systems with long-range interactions in the limit N — > +00 [l2j ■ 

The mean field Kramers equation (|83|) monotonically decreases the free energy functional (|34[) which plays the role 
of a Lyapunov functional. Indeed, a simple calculation gives 

'--/K^fM"**'^ (84) 



For a steady state, F = 0, the term in parenthesis in Eq. (|84|). which is the diffusion current in the mean field Kramers 
equation (|83|) . vanishes. Since d/dt = 0, the advective term (l.h.s.) in Eq. (|83[) must also vanish, independently. From 
these two requirements, we find that the stationary solution of the mean field Kramers equation is the mean-field 
Maxwcll-Boltzmann distribution (f28|) . Therefore, a distribution function /(r, v) is a steady state of the mean field 
Kramers equation iff it is a critical point of free energy (|34p at fixed mass. Furthermore, it can be shown that a steady 
state is linearly dynamically stable with respect to the mean field Kramers equation (|83p iff it is a (local) minimum 
of F at fixed mass. This is consistent with the condition of thermodynamical equilibrium (see Sec. Ill D[) . 



C. Long-range and short-range interactions: Dynamical Density Functional Theory (DDFT) 



In the presence of long and short range interactions, we can simplify the exact Kramers equation (|81 [) by making 
the approximation 



f 2 (r, v, r', v', t) = /(r, v, t)f(v', v', t)p 2 (r, r', t)/p(r, t)p{r\ t). 



This yields 



Of , . df 



dt 



dr p(r,t) 



p 2 (r,r , ,t)S7u{\r-r'\)dr l 



df 



dv dv dv 



df , d fk B _Tdf 



(85) 



(86) 



Using furthermore approximation (|70[) for the short-range interaction and approximation (|66[) for the long-range 
interaction, we get 



df df 6F ex df ^ df d ( k B T df 

ot or dp ov av ov \ m dv 

which is closed. The steady state of this equation is 

/(r,v)=A+(^' + T). 
The integration over the velocity returns the density distribution (|40l) . 



(87) 



V. HYDRODYNAMICS OF BROWNIAN PARTICLES IN INTERACTION 



We now develop a hydrodynamical theory of Brownian particles with long and short range interactions, generalizing 
the theory presented in fill l24l |25| . 
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A. Damped Jeans equations 

Taking the hydrodynamic moments of the exact Kramers equation (|81[) and proceeding as in he mm , we obtain 

^+V-(pu) = 0, (89) 



<9u / ^ 
^ + (u.V)u 



dP f 

~dx^ ~ J P2(r ' r '' *) Vu d r ~ r 'D dr ' ~ P^® ext ~ ( 9 °) 



where p(r,t) = J fdv is the density, u(r, t) = (1/p) J fvdv is the local velocity, w = v — u(r, t) is the relative 
velocity and Pij = J fuiiWjdv is the kinetic pressure tensor. We also recall that the kinetic pressure is defined by 
Pkin(r,t) = ^ J fw 2 dv. For £ = 0, and in the absence of short-range interactions, Eqs. (|89|) - (|90|) reduce to the 
equations obtained by Maxwell in his theory of gases [55, 5j| and by Jeans in the context of stellar dynamics [571 l58j. 
Equations (|59")) - (jM)) will be called the damped Jeans equations (24J. Using the equation of continuity ([BT)|) . we obtain 
the identity 



^ + (u.V)u 



dt {pu) + d^ {pUiU ^ (91) 



On the other hand, proceeding as in [24j |. we obtain the virial theorem 

1 2§ + ¥^- 2e+n+VLR+VsR+v ^ (92) 

where I = J pr 2 dr is the moment of inertia, 9 = | J pu 2 dr is the macroscopic kinetic energy and H = d J pki n dr is 
equal to twice the microscopic kinetic energy (we have K = O + ^11 where K = | J fv 2 drdv is the kinetic energy). 

B. Strong friction limit: Smoluchowski equation 



The exact Smoluchowski equation (|65[) can be derived from the exact Kramers equation (f81j) in the strong friction 
limit £ — > +oo. Considering the r.h.s. of Eq. ([8Tjl. we note that, for £ — > +oo, the velocity distribution is Maxwellian: 



/(r,v,i) = p (r,t)e-^ + OiC 1 )- (93) 

This implies that u = 0(1/0, Pij = (pk B T/m)b V3 + 0(1/0 and Pkin = pk B T/m + 0(1/0- Therefore, to leading 
order in l/£, the damped Jeans equation reduces to 

pu ~ -1 {^^P + J P2(r, r', f)Vu(|r - r'|) dr' + p\7<P ext ^j . (94) 

Inserting Eq. in the continuity equation we obtain the exact Smoluchowski equation (|65[) . This approach 
shows that, for £ — > +oo, the velocity distribution is Maxwellian and the evolution of the spatial density p(r,t) is 
governed by the exact Smoluchowski equation Eq. For £ — > +oo, we have — > and II — > dNk B T, so that the 

virial theorem (JM]) becomes 

= dNk B T + V LR + V SR + Vext- (95) 



This relation can also be obtained directly from the exact Smoluchowski equation (|65|) [241 ] . 

C. Local Thermodynamic Equilibrium (LTE) approximation: damped Euler equation 

The damped Jeans equation (|90l) is not closed since the pressure tensor depends on the next order moment of 
the velocity. Following (m , l24l |25|. we propose to close the hierarchy by making a local thermodynamic equilibrium 
(L.T.E.) approximation: 

Wr.v.t) = ^ d/2 p (r,t)e-^-^ 2 . (96) 
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The distribution function (|96|) minimizes the free energy (|35[) for a given value of density p{r, t) and local velocity 
u(r, t). With the LTE approximation, the pressure tensor takes the form 



Pij = p(r,t) 6i 



(97) 



The kinetic pressure is pk 
du 



P 



Of 



+ (u- V)u 



pfc^T/m. Substituting this result in Eq. (|90[) . we obtain the damped Euler equation 
k B T 



Vp- / p 2 {r,r' ,f)Vu(|r - r'|) dr' - pV<5> ext - &u. 



Using II = dNk B T, the virial theorem (|92[) takes the form 

ld 2 I 1 dl 



2 dt 2 



2^dt = 29 + dNkBT + ^ + ^ + Vext - 



For £ 



-co, we can formally neglect the advective term (l.h.s.) in Eq. (|98[) and we obtain 

'k B T. 



pu 



V/7+ / p 2 {r,r' ,t)Vu(\r - r'\) dr' + pV<P e 



(98) 



(99) 



(100) 



Inserting Eq. (|100[) in the continuity equation (|89[) . we recover the exact Smoluchowski equation (|65[) . However, we 
stress that this cannot be considered as a derivation (even formal) of the exact Smoluchowski equation, unlike the 
derivation of Sec. IV Bl since the damped Euler equation (|9"8")l is heuristic. Indeed, there is no rigorous justification 
of the local thermodynamic equilibrium (L.T.E.) approximation. Accordingly, it does not appear to be possible to 
rigorously derive the damped Euler equation (|9"8"j) from the exact Kramers equation (|5Tj) . 

Remark: The relevance of the LTE approximation has been studied in Ref. 25| in the case of very simple Brow- 
nian models. The conclusion is that the LTE approximation is reasonable and could be improved by replacing the 
temperature of the bath T by a time dependent temperature T(t). 



D. Cattaneo equation 

Using identity (|9"Tj) . the damped Euler equation can be rewritten 
d k T f 

— (pu)+V(pu®u) = — Vp- p 2 {r,r',t)Vu(\r-r'\)dr' - pV^ ext -^pu. (101) 

at m J 



This equation is hyperbolic. If we neglect the inertial term (l.h.s.) in Eq. (|101|) and substitute the resulting expression 
for pu in the continuity equation (|89|). we obtain the exact Smoluchowski equation (|65[) that is parabolic (this is valid 
in a strong friction limit £ — > +oo with the limitation indicated at the end of Se c. IV CI) . The Smoluchowski equation 
neglects memory effects and leads to infinite speed propagation. Following |25|, we can obtain a simplified 
hyperbolic model taking into account memory effects and having a finite speed propagation. Indeed, if we only 
neglect the nonlinear term V(pu <g) u) in Eq. (jlOip . we obtain 



^(pu) = -^lv P - ( p 2 (r, r', t)Wu(\r - r'|) dr' - P V$ ext - frm. (102) 
at m J 



This approximation is exact in the linear regime close to equilibrium where |u| — > [lg. Taking the time derivative 
of Eq. (f8"9")) and substituting Eq. (|102[) in the resulting expression, we find that 



dt 2 ^dt 



Vp + y p 2 (r,r',t)Vu(|r-r'|)dr , + pV* es t) . (103) 



This equation, which is second order in time, is analogous to the Cattaneo equation |59l |. or to the telegraph equation, 
which generalizes the diffusion equation by introducing memory effects and a finite speed propagation. Using — > 0, 
the virial theorem (|9"9"|) takes the form 

1 d 2 I 1 dl 

^^- 2 - + ^=dNk B T + V LR + V SR + V ext . (104) 
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E. Long and short-range interactions 



We now consider the case where the potential of interaction is of the form u = Ulr + usr- We treat the long-range 
interaction in the mean field approximation (|66[) and the short-range interaction with the approximation (|70[) . 
(i) The damped Jeans equations ([89l - (l9T)l) become 



^+V-(pu) = 0, 



(105) 



3u , m 

^ + (u.V)u 



dPu „5F ex „ x „^ 
oxj dp 



The virial theorem is given by Eq. (|92j) where Vlr is replaced by V and Vsr by Vd_ft defined by Eqs. ([32 
(ii) In the strong friction limit £ — > +oo, we obtain the generalized Smoluchowski equation 



(106) 
and (1431). 



<9p 

7# 



v • 



1 /fe B T. 



£ V m 



-Vp + pV- 



.6F„ 



5p 



pV$ + pV$ e 



The virial theorem is given by Eq. (|95[) where Vlr is replaced by V and Vsr by Vdft defined by Eqs. ([32 
If the free energy is of the form (|46|) . then using the relation ([48]) . we get 



(107) 
and (1431). 



9p 

7# 



= v • 



The virial theorem is given by Eq. (|95[) where Vl,r is replaced by F and Vsr by Vdft defined by Eqs. ([32 
It can be written explicitly 

¥f t = d [ pdr+v+v - t - 



(108) 
and (1541). 



(hi) The damped Eulcr equation (|98|) becomes 
<9u 



at 



+ (u- V)u 



k J*v P -pV 5F ' 



5p 



The virial theorem is given by Eq. (|99[) where is replaced by V and Vs^ by Vdft defined by Eqs. 
Note that Eq. (|110[) can be written in terms of the free energy ([35]) as 



du 

~dt 



+ (u • V)u 



T7 SF , 

-pV- £pu. 

op 



This equation satisfies an ii-thcorcm for the total free energy 

F tot [p,u] = J pydr + J p<f> ext dr+^ J p$ dr + k B T J £■ In £dr + F ex [p}. 



For purely long-range interactions, this functional can be deduced from the free energy Q34|) by using 
approximation (|M|) . After some calculations [Hj], we get 



Ft, 



(,pu 2 dr < 0. 



(109) 

(110) 
and gSJ). 

(Ill) 

(112) 
the LTE 

(113) 



For a steady state, F{ j = 0, we obtain u = 0. Then, Eq. (|111[) implies that 5F/5p is uniform leading to Eq. (|37[) . 
Therefore, a steady state of the damped Euler equation is a critical point of free energy at fixed mass. Furthermore, 
it is dynamically stable iff it is a (local) minimum of F at fixed mass M (more precisely, the same results as those 
described at the end of Sec. IIII Cl ean be obtained). If the free energy is of the form (|4l))) . then using the relation 
we can write the damped Euler equation as 



- + (u.V)u 



(114) 
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In that case, the total free energy is 

,2 



F tot [p,u] = J p^-dv + J p<S> ext dr+^ J p$dr + j p J ^J-dp'dr. 



(115) 



The virial theorem is given by Eq. (|99p where Vlr is replaced by V and Vsr by Vdft defined by Eqs. (|32[) and (|54[) . 
It can be written explicitly 

l§ + ¥T t = 2G + d J pdr + v + v ^ ( 116) 

In the strong friction limit £ — > +oo, the preceding equations reduce formally to those obtained in (ii). However, 
as we have already indicated at the end of Sec. IV C[ this is not the correct way to justify these equations since the 
damped Euler equation (|110l) is based on a LTE approximation that has no rigorous foundation, 
(iv) The Cattaneo equation (|103|) becomes 

d 2 p dp fk B T 5F ex \ , s 

It can be written in terms of the free energy ([3"5"|) as 

d 2 p dp ( SF 



at* '*i = v rVJ (118) 

The virial theorem is given by Eq. (|104j) where Vlr is replaced by V and Vsr by Vdft defined by Eqs. (|32f) and 
(143p . If the free energy is of the form (|46|) . then using the relation (j48| . we can write the Cattaneo equation as 



^+^=V-(Vp + pV$ + pV$ ext ). (119) 

The virial theorem is given by Eq. (|104[) where Vlr is replaced by V and Vgfl by Vdft defined by Eqs. (|3^|) and 
(f54|) . It can be written explicitly 

l§ + ¥ft= d J pdr+v+v ^ (i20) 

Remark: starting from Eq. (|87[) . we could also derive a system of hydrodynamic equations for the density p(r,t), 
the local velocity u(r, i) and the kinetic temperature T(r,t). They are given by Eqs. (75)-(77) of with $(r) 
replaced by $(r,t) + $ e *t(r) + f|. 



F. Stochastic kinetic equations 

When there exists metastable states (local minima of free energy), and when the number of particles is not too 
large (77j . it is important to take fluctuations into account. In that case, the preceding deterministic equations must 
be replaced by stochastic equations including a noise term. These equations involve the coarse-grained density p(r, t) 
and coarse-grained distribution function /(r, v, t) that are spatial and/or time coarse-grained averages of the density 
operators pdijc,t) and /d(r,v, t), instead of the average density p(v,t) or average distribution function f(r,v,t) (see 
Appendix |B|) . For finite iV systems, or when we are close to a critical point, the fluctuations can induce random 
transitions from one metastable state to the other. The system will visit these minima randomly and pass from one 
state to the other (of course, the global minimum of free energy is the most frequently visited). On the other hand, 
for large N, or when we are far from a critical point, the metastable states have extremely long lifetimes (larger in 
practice than the duration of the experiment), scaling like e N (49[, and they are as much relevant as fully stable states. 
In that case, we can use the deterministic equations of the previous sections. The relation between deterministic and 
stochastic equations, and the importance of metastable states, are further discussed in Appendix iBl 



The stochastic coarse-grained Kramers equation is 16j : 
I + " • I " I '/ ^ r', V, f )V„(|r - - ^ • + 7v) + I ■ (v^7Q(r, v, t ) 



(121) 
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where Q(r, v, f) is a Gaussian white noise satisfying (Q(r, v, £)) = and (Q(r, v, i)Q(r', v', t')) = S(r — r')6(v — 
v')S(t — t'). We can now generalize the derivation of the hydrodynamical equations given in Ref. [Tr| . 
The stochastic damped Jeans equations are 

g + V-(pu) = 0, (122) 



du 

~dt 



+ (u • V)u 



dP 

dx 



^ - Jp^(r, r 7 , *)Vu(|r - r'|) dv' - pV$ eTt - £pu - y/2^k B TpR(r, t). 



(123) 



Using approximations similar to those made in Sees. IIIIBI and 1111 CI but now applying to the coarse-grained distri- 
bution functions, we replace the integral involving the correlation function by 



SF — 
pl(r, r', t)Vu(r, r') dv' = p(r, t)V— ^[p(r, t)] + p(r, t)V$(r, t), 

dp 



(124) 



where Fezlp] is the equilibrium free energy functional and $(r, t) the mean field potential (|68p determined by the 
coarse-grained density p(r, i). This leads to the stochastic damped Jeans equation 



^ + (n.v)n 



dPj. 



.SF, 



- PV— ^ - pV$ - pV$ e:ct - £p u - v /2^fc s TpR(r, i). 
dp 



In the strong friction limit, we get the stochastic coarse-grained Smoluchowski equation 

V ■ 



&p 



^ = V- -pV— 



0f 



1 (5F 



l 2k B Tp 



R 



(125) 



(126) 



where F[p] is the free energy functional (|35|). 

If we make the LTE approximation, we get the stochastic coarse-grained damped Eulcr equation 



du 

~dt 



(u • V)u 



PV— - £pu - ^2£fc B TpR(r, t). 
dp 



Finally, the stochastic coarse-grained Cattanco equation is 



The virial theorem associated with the stochastic damped Eulcr equation (|127j) is 



(127) 



(128) 



(129) 



where r)(t) is a Gaussian white noise satisfying (r](t)) = and (r](t)r](t')) = S(t — if). Using Eqs. ([33)) , ([3"2l and 
the second term in the r.h.s. can be decomposed into 



pr ■ dr = d/Vfc B T + V DFT + V + V ext . 
Sp 



The virial theorem associated with the stochastic Smoluchowski equation (|1 26|) is 

^§ = - J pr • dr - y/2tk B TI(t)r,(t), 

and the virial theorem associated with the stochastic Cattaneo equation (|128[) is 

ld 2 I ldl f SF , r— 

2d^ + ^Tt=-J~ pr - v ^ dr -^ m ^- 



(130) 



(131) 



(132) 
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VI. NONLINEAR MEAN FIELD FOKKER-PLANCK EQUATIONS AND GENERALIZED 

THERMODYNAMICS 



Hydrodynamic equations similar to those derived previously, including a long-range potential of interaction and a 
generically nonlinear equation of state p(p) taking into account small-scale constraints, had previously been derived 
[44j in the context of nonlinear mean field Fokker-Planck equations (NFP) and generalized thermodynamics pioneered 
by Tsallis and Plastino & Plastino [4l| (see [38|, [33, |42j for reviews). However, the origin of the nonlinear 
equation of state p(p) is physically different in the two approaches. In the context of generalized thermodynamics, 
the equation of state arises from the non-Boltzmannian nature of the distribution function (and p = puin represents 
the kinetic pressure) while in the theory of fluids developed in Sec. [Vj the distribution function is Boltzmannian and 
the equation of state arises from the two-body correlation function induced by the short-range potential of interaction 
(andp represents the thermodynamic pressure). Despite this fundamental difference, the hydrodynamic (macroscopic) 
equations have the same mathematical form! It is therefore interesting to compare these two approaches in detail in 
order to stress their analogies and differences. 



A. Overdamped model 

Let us consider a system of Brownian particles in interaction in the overdamped limit. We assume that the particles 
interact via a mean field potential $(r, f) given by Eq. (j68]) and that they are submitted to an external potential 
<5> ext {r). For the moment, we ignore small-scale constraints. In that case, the motion of a particle is described by the 
stochastic Langevin equation: 



(133) 



where £ and D arc the coefficients of friction and diffusion and R(t) is a Gaussian white noise. The temperature 
is given by the Einstein relation D = £T (we take the Boltzmann constant and the mass of the particles equal to 
unity). The corresponding Fokker-Planck equation is the mean field Smoluchowski equation (|67p. In order to take 
into account microscopic constraints that affect the motion of the particles, Kaniadakis [43| has proposed to modify 
the form of the transition probability from one state to another. This kinetical interaction principle (KIP) can take 
into account exclusion or inclusion constraints that enhance or inhibit the transition. This can model for example 
quantum effects, close packing effects, steric hindrance... Let us define the transition probability of a particle from 
position r to position r' by (see [43| and Sec. 2.11 of [39] for details): 

7r(r -> r') = w(r, r' - r)a[p(v, t)]b[p(v\ t)], (134) 

where w(r,r' — r) is the transition rate that only depends on the nature of the interaction between the test particle 
and the bath, and a(p) and b{p) are positive functions. Linear kinetics corresponds to a(p) = p and b(p) = 1: the 
transition probability is proportional to the density of the starting state and independent on the density of the arrival 
state. It leads to the ordinary Fokker-Planck equation (|67[) . Here, we assume a more general dependence on the 
occupancy in the starting and arrival states. This creates a biais with respect to the ordinary situation. Using a first 
neighbor ap pro ximation and an extension of the Kramers- Moyal expansion based on the transition probability f|L 34[) , 
Kaniadakis [43j obtains a nonlinear Fokker-Planck equation of the form 

= V ■ [Th{p)V P + fl(p)V($ + , (135) 

where the functions h(p) and g(p) are related to the bias a(p) and b(p) in the transition probabilities by 

g(p) = a(p)b(p), h(p)^b(p)a'( P )-a(p)b / (p). (136) 
The generalized free energy associated with the NFP equation (|135p is 

F[p] =E-TS=^J pt>dr + J p<P ext dr + T J C{p) dr, (137) 



where 



C(p)dv, C "(p) = ^-, (138) 
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is a "generalized entropy" determined by the ratio of the functions h(p) and g(p) [38|, [39j, |43j, |44|, [6(J . We shall assume 
that C is convex (i.e. C" > 0). In the absence of microscopic constraint, h(p) = 1 and g(p) = p. In that case, Eq. 
(|135| reduces to the mean field Smoluchowski equation (|67|) . the entropy (|138p reduces to the Boltzmann entropy 
S = — J plnpdr and the free ener gy (1 1371) reduces to the Boltzmann free energy (f3"3"|) . We note that the NFP equation 
(|135p can be written in the form 38|, 39j ] : 



dp 
dt 



-^ p)V -p 



This equation satisfies an ii-theorem: 



F 



g(p) v 



.SF 



Sp 



dr < 0. 



(139) 



(140) 



For a steady state, F = 0, Eq. (|140[) implies that SF/Sp is uniform. Therefore, a steady state of the NFP equation 
(11391) is a critical point of free energy at fixed mass. Furthermore, it is dynamically stable iff it is a (local) minimum 
of F at fixed mass M (more precisely, the same results as those described at the end of Sec. IIII Cl ean be obtained). 
Writing 5F + aTSM = 0, where a is a Lagrange multiplier, we find that the steady states of the generalized mean 
field Smoluchowski equation (|1 35[) are given by 



p(r) = (CT 1 + *«xt) -a] ■ 

In the following, we shall assume that g(p) = p. In that case, the NFP equation (|135[) becomes 

^ = V • [TpC"( P )Vp + P V($ + $ ext )] . 
It can be written in the form of a generalized mean field Smoluchowski equation 

f^ = V-[Vp + pV(* + $ext)], 



with a barotropic equation of state p(p) given by 



V '(p)=TpC"{p). 



(141) 



(142) 



(143) 



(144) 



Equation (|143p is mathematically equivalent to Eq. (|72p. Its steady state is given by the condition of hydrostatic 
equilibrium (|52[) . Since C(p) is convex, we find that p'{p) > 0. A first integration gives 



C( P ) 



p(p) = Tp 2 
A second integration leads to the identity 

TC( P ) = p 

Therefore, the free energy (|137p can be rewritten 



= T[C'(p)p-C{p)\. 



P P(P1), 
— o— dpi. 
Pi 



F[p] = ^Jp<S>dr + J p^ ext dr + J pj 



— —dpidv. 



Pi 



(145) 



(146) 



(147) 



It coincides with the free energy (|46p describing a fluid with weak density gradients. Therefore, the two approaches 
(DDFT and generalized thermodynamics) lead to similar equations but for different reasons. We finally note that the 
NFP equation (|143[) can be derived from the generalized stochastic process 



dv 




$p(pM)) 



R(t), 



(148) 



in which the noise term explicitly depends on the density of particles around the particle under consideration. This 
is a phenomenological manner to take into account microscopic constraints that can affect the motion of the particles 

SIM El. 
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B. Inertial model 



We now extend the previous approach in phase space in order to take into account the inertia of the particles. In 
the absence of microscopic constraints, the mean field dynamics of Brownian particles is described by the Langevin 
equations 

§ = v. (149) 

^ = -f v - V$- V$ ext + V2DR(t). (150) 

The ordinary Fokker-Planck equation associated with these equations is the mean field Kramers equation (|53")) . Mod- 
eling microscopic constraints with the KIP J43|, and defining the transition probability by 

7r(r, v -> v') = «;(r, v, v' - v)o[/(r, v, t)]b[f(v, v', t)], (151) 

Kaniadakis [43| obtains a nonlinear Fokker-Planck equation of the form 



dt dr ext dv dv 



Th(f)^+g(f)v 



(152) 



where the functions h(f) and <?(/) are related to the bias a(/) and b(f) in the transition probabilities by 

g(f) = a(f)b(f), h(f) = b(f)a'(f) - a(f)b'(f). (153) 
The generalized free energy associated with the NFP equation (|152[) is 

F[f] =E-TS=^J fv 2 drdv +lj p ^dr + J p$ ext dr + T J C(f) drdv, (154) 



C(/)drdv, C"(/) = ^, (155) 

9(f) 



where 



is a "generalized entropy" determined by the ratio of the functions h{f) and g(f) |38j, |39|, |43j, |4J, |60( • In the absence 
of microscopic constraints, a(f) = f and b(f) = 1, implying /i(/) = 1 and g(/) = /. In that case, Eq. ()152[) reduces 
to the mean field Kramers equation ((53")) , the entropy (|155[) reduces to the Boltzmann entropy S = — J f In / drdv 
and the free energy (11541) r educes to the Boltzmann free energy (p4|) . We note that the NFP equation (|152[) can be 
written in the form [38l . |39( : 



9 ^ N 2 



This equation satisfies an _ff-thcorcm: 

F=-Jt9(f)^f I ) drdv<0. (157) 

For a steady state, F = 0, Eq. (|157|) implies that SF/Sf is independent on v so that the current in the NFP equation 
(I15GI) vanishes. Since dtf = 0, the inertial term in Eq. (|156j) must also vanish, independently. From these two 
requirements, we find that 5F/8f is independent on r. As a result, 6F/6f is constant. Therefore, a steady state of 
the NFP equation (| 1 56(1 is a critical point of free energy at fixed mass. Furthermore, it is dynamically stable iff it is 
a (local) minimum of F at fixed mass M (more precisely, the same results as those described at the end of Sec. IIII CI 
can be obtained). Writing SF + aTSM = 0, where a is a Lagrange multiplier, we find that the steady states of the 
generalized mean field Kramers equation (|152[) are given by 



/(r,v) = (C")-M-/3 



— + $(r) + $ext(r) 



(158) 
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Introducing the density p = J f dv and the kinetic pressure P = 2 f f v2 dv, and using Eq. (|158p . we find that 
p = p{f3® tot (r) + a] and p = p[/3$ tot (r) + a], where $ tot = $ + $ ext . Eliminating /3$ tot (r) + a between these 
two expressions, we find that the equation of state at equilibrium is barotropic: p = p(p). We emphasize that the 
pressure here defined is the kinetic pressure and that the equation of state p = p(p) is completely determined by 
the function C(f), hence by the bias a(f) and b(f) in the transition probabilities. In the usual (Boltzmann) case 
where C(f) = / In/, we get the isothermal equation of state p = pT (linear) but for more general functions C(f), the 
equation of state p{p) is nonlinear. Finally, it is easy to check (see, e.g., [33]) that Eq. (|158p implies the condition of 
hydrostatic equilibrium 



Vp + pV$ + pV$ e 



0. 



(159) 



We therefore obtain the same result as in Eq. (j52"]) but for a fundamentally different reason. In Eq. (f5"2"j) . the quantity 
p designates the thermodynamical pressure p = pid +p e x, where pid = pfcsT/m is the ideal pressure (coinciding with 
the kinetic pressure) and p ex is the excess pressure taking into account short-range interactions. In that approach, the 
velocity distribution is Maxwcllian leading to the isothermal gas law pid = phsT/m and the excess pressure comes 
from correlations. On the other hand, in Eq. (|159|) . the quantity p designates the kinetic pressure. In that approach, 
correlations are ignored (or taken into account implicitly in the KIP) but the velocity distribution is non-Maxwellian 
leading to a nonlinear equation of state p(p). 

In the following, we shall assume that g(f) = f. In that case, the NFP equation (|152p becomes 



df 
dt 



or 



df _ d 



T/C"(/)^ + /v 



(160) 



It can be derived from the generalized stochastic process 



dv 

dt 



(161) 



^ = -£v - V$ - V<t> ext + J 2D J 



C(f) 



f 



R(t), 



(162) 



where the noise explicitly depends on the distribution of particles around the particle under consideration (in phase 
space) HIM HI- 



C. The strong friction limit 



In order to stress the differences with the results of Sec. IIV1 we shall recall the derivation of the generalized 
Smoluchowski equation (|143[) from the generalized Kramers equation (|160[) in the strong friction limit £ — > +00 [39( |. 
The first two moments of the hierarchy of hydrodynamic equations associated with Eq. (|160|) are 



dp 
dt 



V ■ (pu) = 0. 



(163) 



- + (u.V)u 



dPu 

_ pV $ - pV$ ext - £pu, 

dxj 



(164) 



where the quantities have the same meaning as in Sec. IV Al We now consider the strong friction limit £ — > +00 with 
fixed T . Since the term in parenthesis in Eq. (|160[) must vanish at leading order, wc find that the out-of-equilibrium 
distribution function /o(r, v, f) is given by 



/o(r,v,t) = (C")- 1 |-/3 y+A(r,t) j+O^ 1 ), 
where A(r,i) is a constant of integration that is determined by the density according to 

p(r,t) = [ f dv = p[X(r,t)}. 



(165) 



(166) 
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Note that the distribution function fo is isotropic so that the velocity u(r, t) = 0(£ 1 ) and the pressure tensor 
Pij = pSij + 0(£ _1 ) where p is given by 



f v 2 dv = p[X(r,t)}. 



(167) 



Eliminating X(r,t) between the two expressions (|166|) and (|167l) . we find that the fluid is barotropic with an equation 
of state p = p(p), the same as in equilibrium (see Sec. IVI Bj) . Now, considering the momentum equation (|164j) in the 
limit £ — > +00, we find that 



(168) 



Inserting this relation in the continuity equation f|163[) , we obtain the generalized mean field Smoluchowski equation 

1 



dp 

m 



(Vp + pV$ + pV$ ea;t ) 



The free energy associated with this equation is 



2 dpidr + — I p$dr 
Pi 2 



p^extdr. 



(169) 



(170) 



It can be deduced from the free energy (| 1 54(1 by using Eq. (|165|) to express F[f] as a functional F[p] = F[fo] of the 
density (see [39| for the details of calculation). This leads to the same equations as in Sec. IVBI but. as explained at 
the end of Sec. IVIB1 the reason is fundamentally different. 

Remark: the generalized Smoluchowski equation can also be derived from the generalized Kramers equation by 
using a Chapman-Enskog expansion (62| . In that case, it is possible to consider generalized Kramers equations of 
the form (|152p with arbitrary g(f). This leads to generalized Smoluchowski equations of the form (|169p where £ now 
depends on position and time. 



D. Damped Euler equation 



We can also derive a damped Euler equation similar to the one obtained in Sec. IV CI To that purpose, we close 
the damped Jeans equation (|164|) by using a LTE approximation fjij : 



/ iTB (r ! v ! i) = (C")- 1 |-^ 
where A(r,i) is determined by the density according to 



(v-u(r,£)) 2 



+ A(r,t) 



P(r, t) = J fLTEdv = p[\(T, t)] 



(171) 



(172) 



The distribution function (|171[) minimizes the free energy (|154[) for a given value of the density p{r,t) and local 
velocity u(r, t). With the LTE approximation, the pressure tensor takes the form = pSij where p is given by 



p(r,t) = i J / LT B(v-u(r,i)) 2 dv=p[A(r,t)]. 



(173) 



Eliminating A(r, t) between the two expressions (|1T2|) and f|l T3[) . we find that the fluid is barotropic with an equation 
of state p = p(p), the same as in the preceding sections. Substituting these results in Eqs. (|163p and (|164j) . we obtain 
the damped Euler equations 



| + V.(pu)=0. 



(174) 



- + (u.V)u 



(175) 
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The free energy associated with these equations is 

F tot [p,u} = J p\dv + J p J f ' ?kftd Pl dr+^ J p$dr + J p$ ext dr. (176) 

It can be deduced from the free energy (|154|) by using Eq. (|171j) to express F[f] as a functional F[p, u] = FI/lte] of 
the density and local velocity. This leads to the same equations as in Sec. IV CI but, as explained at the end of Sec. 
I VI Bl the reason is fundamentally different. 



E. Generalized Cahn-Hilliard equations 



Let us assume that the long-range potential u(\r — r'|) is screened on a distance that is large with respect to the 
microscopic length but short with respect to the system size. Therefore, we assume that the generalized mean field 
Smoluchowski equation (|142|) remains valid, but that we can simplify the potential &(r,t) given by Eq. (|68|) . Setting 
q = r' r and writing 



$(r,t) = J u(q)p(r + q,t)dq, 
we Taylor expand p(r + q, t) up to second order in q: 

dp 1 v-^ d 2 p 
p(r + q, i) = p(r, t) + J2 a£<R + ^ E 

i 

Substituting this expansion in Eq. (|177[) . we obtain 



'■j 



dxidx 



*(r,t) 



-ap(r,t) - -Ap(r,t), 



(177) 



(178) 



(179) 



with a = —Sd / +O ° u(q)q d ~ 1 dq and b — Jo + °° u {<l)<l d+1 dq- Note that I = (b/a) 1/>2 has the dimension of a length 

corresponding to the range of the interaction. Substituting Eq. (|179[) in Eq. (|137|) , we can put the free energy in the 
form 



F[p] 



(Vp) 2 + V(p) 



dr. 



where V is the effective potential 



V(p) 



IT 



C(p). 



In that case, Eq. (|139[) can be rewritten 



^ = -V-[Ag(p)\7(Ap-V'(p))} 

with A = 6/(2£). The steady state of Eq. (| 139[) or (|182l) corresponds to a uniform SF/Sp = —aT leading to 

2aT 



Ap = V'(p) + 



(180) 



(181) 



(182) 



(183) 



Equation (|182[) share some analogies (but also crucial differences) with the Cahn-Hilliard equations [63( ; see discussion 
in (l6| . Coincidentally, the case g(p) = 1 and C(p) = p 4 gives an equation that is formally equivalent to the Cahn- 
Hilliard equation with V(p) = ^&{-^ — p 2 ) 2 ■ On the other hand, in the classical (Boltzmann) case g(p) = p and 
C = p In p, Eq. (|182l) takes the form 



(T-ap)Vp--pV{Ap) 



(184) 



involving a density dependent diffusion coefficient D{p) = |(T — ap) and an effective potential $ e // = — j^^P- 

Remark: as discussed at the end of Sec. Ill E 51 this gradient expansion can also be performed in the case of systems 
with short-range interactions described by the free energy functional (|56|) . provided that we use the correspondence 
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VII. CONCLUSION 



In this paper, we have developed a general kinetic theory of Brownian particles with long and short range in- 
teractions. To close the BBGKY-like hierarchy, we have used the mean field approximation (|5t?|) for the long-range 
interaction Q and the standard approximation (|70[) used in the theory of fluids for the short-range interaction [35| . We 
have included these approximations in the general kinetic and hydrodynamic equations of Brownian particles derived 
in [3 EB, HE] ■ in the overdamped limit, this leads to the generalized mean field Smoluchowski equation (fT2"j) which 
includes a mean field potential due to long-range interactions and a generically nonlinear pressure due to short-range 
interactions. More general equations taking into account inertial effects have also been obtained. Therefore, our 
kinetic theory justifies from a microscopic model the basic equations that have been introduced phenomenologi cally 
to describe various systems with long and short range interactions such as self-gravitating Brownian particles [261 ] . 




chemotaxis of bacterial populations |2jj and colloidal particles with capillary interactions |28|. The justification of 
these models from a kinetic theory starting from microscopic processes was the main goal of this paper. 

We have also found that the same hydrodynamic (macroscopic) equations are obtained from nonlinear mean field 
Fokkcr-Planck equations based on generalized thermodynamics. In both cases, the nonlinear pressure takes into 
account microscopic constraints that affect the dynamics of the particles at small scales. However, the origin of this 
pressure is different. In the classical theory of fluids [HI, [Hj], the distribution function is Boltzmannian and the 
nonlinear pressure is due to the two-body correlation function induced by the small-scale potential usr. Different 
methods have been developed in the theory of fluids to obtain the expression of the pressure law p = p(p 7 T) depending 
on the short-range interactions. In the generalized thermodynamics approach [HI, HH, |42| , the nonlinear pressure arises 
from the non-Boltzmannian nature of the distribution function due to the bias in the transition probabilities from 
one state to the other [43| . It is interesting to observe that the hydrodynamic (macroscopic) equations coincide while 
the microscopic models are fundamentally different. 



Let us consider a spatially homogeneous fluid enclosed within a container of volume V. The pressure p = p(p, T) 
is a function of the density p = M/V and temperature T. We assume that the fluid is in contact with a thermal 
bath imposing the temperature. Since T is fixed, the pressure is barotropic so that p = p(p). Introducing the 
free energy F = E — TS and using the first law of thermodynamics dE = —pdV + TdS, we obtain the identity 
dF = —pdV where we have used dT = 0. This can be rewritten dF = —pMd(l/p) = (p/ p 2 )Mdp and, in integral 
form, F = M J P {p/p 2 ) dp. Introducing the free energy per unit volume / = F/V, we obtain f{p) = p J P {p(p)/p 2 ) dp. 
This relation remains valid locally in an inhomogeneous fluid provided that there are no strong gradients of density 
(e.g., the fluid is sufficiently far away from an interface). This leads to the following expression of the free energy 
functional 



In this Appendix, we discuss the importance of metastable states in the kinetic theory of systems with long (and 
short) range interactions. This will help us to better understand the relation between deterministic and stochastic 
kinetic equations. 

As emphasized by Archer & Rauscher [36J, we must distinguish three types of density fields: (i) the density 
operator pd(r,t) = J^iLi m6(ri(t) — r) which is made of a sum of Dirac peaks coinciding with the exact positions of 
the particles, (ii) the ensemble average density field p(r, t) = (Xa=i n}S(ri(t) — r)) = NmPi(r, t), and (iii) the coarse- 
grained density field p(r, t) which can be viewed either as a spatial [64 | or time [36j average of the density operator 
Pd(r, t) — m ^( r i(t) ~ r). Note that this distinction is standard in equilibrium statistical mechanics. Historically, 

it first appeared in Boltzmann's combinatorial analysis. A microstate is characterized by the density pd(r, t) specifying 
the exact position of all the particles while a macrostate is characterized by the coarse-grained density p(r) such that 
~p(r)dr gives the number of particles in a macrocell [x,x + dx] x [y,y + dy] x [z,z + dz] irrespectively of the exact 
positions of the particles in the cell. At statistical equilibrium, the density probability of the coarse-grained density 
field ~p(r) is 




Appendix A: Expression of the free energy in the absence of strong gradients 




(Al) 



Appendix B: Some comments about the importance of metastable states 



1 



e - pF[p] 8(M - M[p}), 



(Bl) 
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where F[p] is the free energy (|35|) . The normalization condition J P eq [p] Up — 1 leads to the expression of the partition 
function Z(fi) = J e~^ F ^S(M - M\p])Trp. The ensemble average density /O eg (r) corresponds to the average value 
of ~p(r), i.e. Peq(r) = / Peq\fi\pT^P- At the thermodynamic limit, it coincides with the most probable value of ~p{r). 
Therefore, p eq (v) corresponds to the global minimum p g i bai{^) of F[p]. For systems with long-range interactions, this 
has been proven rigorously in [l8| . 

When the free energy F[p] has a unique (global) minimum, the situation is simple. The density probability P e o[p] 
obtained at a given time t from an ensemble of experiments, or the density probability P eq [p] obtained from the time 
series of a unique experiment, coincide and are given by Eq. (|B1|) . On the other hand, the average (or most probable) 
value of p(r) is given by p e q(r) which is the (global) minimum of F[p]. Let us now consider the more complicated 
situation where F[p] has several minima (metastable states). This occurs in particular for systems with long-range 
interactions, like self-gravitating systems [HI EH, and we shall focus on these systems in the following discussion. 
If we fix the time t large enough [78| and consider an ensemble of experiments, we will measure a coarse-grained 
density field p(r) that fluctuates from experiment to experiment. Its density distribution will be given by Eq. (IB It . 
The system will be found most of the time in a minimum of F[p], the global minimum being the most frequent one. 
Recalling that for systems with long-range interactions the free energy is extensive, the equilibrium density probability 
of the coarse-grained distribution (|B1[) can be rewritten 

Pe " m = zJp) e ~ mm5{M ~ M[n) > (B2) 

where f[p] = F[p]/N is independent on N. For N — > +00, the distribution is strongly peaked around the global 
minimum of F[p] at fixed mass, so that an overwhelming majority of configurations with ~p(r) ~ p g iobai(r) we will 
observed. Accordingly, the partition function is dominated by the contribution of the global minimum and we can 
make the approximation Z(j3) ~ e~ pNfm ^^ so that F{fi) = -jhiZ(p) ~ Nf min ((3) where f mm ((3) = f[p 3 iobal}- 
Equivalently, we have P e g[7>] — S(~p — Pgiobal)- However, focussing exclusively on the distribution P eq [p] and on the 
ensemble average p eq (r) may hide the importance of metastable states in the dynamics (79j . To see that, let us now 
consider one experiment and follow the system in time. If N is not too large and/or if we wait long enough, we will sec 
that the system undergoes random transitions from one metastable state to the other. Of course, the global minimum 
of free energy is the most frequently visited. The residence time in a metastable state depends on the barrier of free 
energy with the other minima and is given by the Kramers formula tuf e ~ e^' AF '. If we compute the pdf of ~p~{r,t) 
on a time series over sufficiently long times, we will obtain the distribution (|B1|) . However, we insist on the fact that 
the time on which we calculate the pdf must be extremely long (recall that equilibrium statistical mechanics assume 
ideally that t — > +00). Indeed, since the free energy is extensive, the barrier of free energy scales linearly with the 
number of particles (|AF| ~ N) and the lifetime of a metastable state scales like tuj e ~ e N [49[. Therefore, when N 
is large (it has not to be very large since the dependence of tuf e with ./V is exponential), the metastable states will 
have tremendously long lifetimes! In practice, the system may remain blocked in a metastable state p(r, t) ~ p me ta(r) 
for all the duration of the physical experiment. In that case, we will measure P p hys[p] — S(~p — p me ta) instead of P eq [p]. 
This leads to an apparently non crgodic behavior, although crgodicity holds provided that we wait long enough. 
These arguments show that the limits N — > +00 and t — > +00 do not commute. If we fix ./V and make an experiment 
on a duration t — > +00, we will obtain -P e9 [/5] given by Eq. (|B1[) and p(r) = p g iobai(j)- Alternatively, if we fix an 
interval of time [0,T] that is large but finite and let N — > +00, the system may remain blocked in a metastable 
state during all the duration of the experiment (since its lifetime diverges like e N — > +00). In that case, we will find 
Pphys[p] — S(~p — pmeta) and p(v) = p m eta{^) ■ Having realized that, we can now better understand the relation between 
the deterministic and the stochastic kinetic equations. 

The time evolution of the density operator pd(r,t) = Yli=i TnS(ri(t) — r) is given by the exact stochastic kinetic 
equation (65j : 



dt \£ S Pd 



where 



F d [Pd] = \ [ p d (v,t)u(\r-r'\)p d (v',t)drdr'+ f Pd <5> ext dv + T f ^In^dvdr', (B4) 
2 J J J m m 

is the exact free energy This equation bears exactly the same information as the iV-body dynamics (|59|) and, as 
such, is not very useful for practical applications. The stochastic kinetic equation (|B3[) can be viewed as a Langevin 
equation with a multiplicative noise that vanishes when p d = 0. This expresses the fact that the density cannot 
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fluctuate in regions devoid of particles. The corresponding Fokkcr-Planck equation for the density probability P[pd, t] 
of the distribution pd(r,t) is given by }o6j : 



dP[ Pd ,t] 
dt 



Spd{r,t) 



V -p d {r,t)V 



T 



SF d 



Sp d (r,t) Sp d (r,t) 



P[pd,t] \ dr 



(B5) 



The steady state of this Fokkcr-Planck equation is P eq [pd] oc e l 3Fd ^P d ]S(J Pddv — N) [(55[ which is equivalent to 
the canonical N body distribution (J7J). The ensemble average density p(r,t) = (J2iLi mS(ri(t) — r)}, which is a 
deterministic field, satisfies the exact Smoluchowski equation (|65[) . It can be obtained by averaging the exact stochastic 
equation for pd(r,t) [l|| [H[ or by writing the first exact equation of the BBGKY hierarchy [12, However, this 
equation is not closed and some approximations, whose validity will be discussed below, must be introduced. Finally, 
the coarse-grained density ~p(r, t) is a fluctuating field whose evolution is governed by the stochastic kinetic equation 
(If 2GI) with the free energy (|35[) . This equation can be obtained in a strong friction limit of fluctuating hydrodynamic 
equations [68[ , by coarse-graining the exact stochastic equation (|B3[) [36j or by using the general theory of fluctuations 
of Landau & Lifshitz [(39[ (see Appendix B of [Hj]). The stochastic kinetic equation (|1 26[) can be viewed as a Langevin 
equation for the coarse-grained density ~p(r,t). The density probability P[p,t] of the coarse-grained density ~p{v,t) is 
given by the Fokker-Planck equation [64l . l6q : 



dt 



Sp(r,t) 



V-p(r,t)V 



T 



5F 



5p{r,t) Sp(r,t) 



P[p,t] \ dr 



(B6) 



where F[p] is the free energy (|35[) . The steady state of this Fokkcr-Planck equation is the statistical equilibrium 
state (|Bip foil l68j . The convergence towards this equilibrium state is guaranteed by an if-theorem [7(3| . The density 
p(r,t) is the average value of p(r, t), i.e. p{r 7 t) = J P[p,t]pT>p. In the thermodynamic limit, it coincides with the 
most probable value of p(r,i). As we have previously indicated, the density p(r,t) satisfies the exact Smoluchowski 
equation (|65j) . but this equation is not closed. The stochastic kinetic equation (|126|) . which is closed (unlike the exact 
Smoluchowski equation ([55)) ) and which describes the evolution of a smooth field (unlike the exact stochastic equation 
(|F33|) which describes the evolution of an operator made of Dirac peaks) is essentially "exact" and represents the most 
important equation of the list. 

Let us now discuss the validity of the approximate Smoluchowski equation (|73p in the light of the previous con- 
siderations. If the functional F[p] has a unique (global) minimum, then the solution of the deterministic equation 
(|73p converges towards this minimum for t — > +oo. Therefore, the density p(r,t) tends to the equilibrium density 
Peq(j) and Eq. ([73]) certainly provides a good description of the average dynamics. However, when the free energy 
functional F[p] possesses several local minima (metastable states), the situation is more complicated. It that case, the 
deterministic kinetic equation (|73[) will converge for t — > +oo to one of these minima (local or global), whose selection 
will depend on a notion of basin of attraction. Since Eq. (|73[) is a deterministic equation, the system will remain in 
that state p m eta{^,t) for ever, even if this is not the global minimum of F[p]. Therefore, different initial conditions 
(belonging to different basins of attraction) will lead to different density profiles p(r) for t —> +oo in contradiction 
with the fact that p(r,t) should tend to a unique profile Peqi?) which is the average (or most probable) value of p(r) 
according to the distribution ([Bip . This indicates that the approximate deterministic equation (|73|) is not correct 
when there exists metastable states since it may not converge towards the correct equilibrium state. By contrast, 
the average (or most probable) density profile p(r, t) determined from the solution of the stochastic equation (|126p 
correctly tends towards the equilibrium profile p eq {r) for t —> +oo. This shows that the stochastic equation (If 261) is 
superior to the deterministic equation ([73[) . 

However, for systems with long-range interactions, it has been proven rigorously that the mean field approximation 
is exact for N — > +oo, so that the mean field Smoluchowski equation (|57| should be exact in that limit even if there 
exists metastable states. How can we solve this apparent paradox? The solution comes from the non-commutation of 
the limits t — > +oo and N — > +oo that we have previously indicated [80| . The validity of the mean field Smoluchowski 
equation (p7|) assumes that we fix the interval of time [0,T] (any), then let N —> +oo. In that case, we have seen 
that the lifetime of a metastable state tends to e N — > +oo so that it is larger than the physical interval [0,T]. 
This is the reason why the solution of the mean field Smoluchowski equation can gcnerically converge towards any 
minimum (local or global) of F[p] and stay there permanently (it does not make any difference between fully stable or 
metastable states). Alternatively, if we fix N (any), then let t — > +oo, the evolution of p(r, t) is not described by the 
mean field Smoluchowski equation. In that case, the coarse-grained density p(r, t) undergoes random changes from 
one metastable state to the other and we must consider the stochastic equation (|126p if we want to take into account 
these random changes properly. 

In practice, if N is large, the system may be trapped for all physically relevant times in a metastable state which 
is not the global free energy minimum. Still, this metastable state is fully relevant on a physical point of view. It has 
an extremely long lifetime, scaling like e N [49[, so that it is as much relevant as the fully stable state when N — » +oo. 
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This is a situation of physical ergodicity breaking (although, strictly speaking, the system is ergodic if we wait long 
enough). In that case, the system will not jump to another metastable state in the duration of the experiment, so 
that we can use the mean field Smoluchowski equation (|67j) which is deterministic. On the other hand, if N is not 
too large, or if we are close to a critical point so that the barrier of free energy \AF\ is small we will observe 
random changes from one metastable state to the other in the duration of the experiment. In that case, we must use 
the stochastic kinetic equation (|126j) . Typically, we must use the stochastic kinetic equation (|126[) when the physical 
timcscalc of the experiment is larger that e N or more precisely e N ^ A ^. Therefore, the domain of validity of the mean 
field Smoluchowski equation (p7|) for systems with long-range interactions is clearly established. By contrast, the 
domain of validity of the approximate equation (|73|) for systems with short range interactions remains less clear when 
there exists metastable states since there is no small parameter (like 1/N in the previous case) in the limit of which 
this equation can be rigorously justified. Nevertheless, by analogy, we can argue that this equation can be employed 
when the timescale of the experiment is much smaller than the typical lifetime of a metastable state. 
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